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We present an approach based on a non-Hermitian Hamiltonian to describe the process 
of measurement by tunneling of a phase qubit state. We derive simple analytical expres- 
sions which describe the dynamics of measurement, and compare our results with those 
experimentally available. 
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1. Introduction 

The superconducting phase qubit has significant potential applications in quantum 
information processing due to its convenient functional operations, its capability 
to reduce the effects of n oise and decoherence, and its future utilization in scal- 
able quantum computers The improvement of the fidelity of the qubit state 
measurement still remains one of the main issues. 

In this paper, we propose an approach based on non-Hermitian Hamiltonians 
to describe t lie dynamics of measurement of the superconducting phase qubit state 
| 1 | 2 | 3 | 4 | 5 [617J We obtain simple analytical expressions for the probability of tun- 
neling to the continuum. We com pare our theoretical results with experiments on 
superconducting phase qubits and demonstrate that our theoretical predictions 
are in a good agreement with the available experimental data. 



I 



I 



18, 2010 17:39 WSPC/INSTRUCTION FILE IJQFPhQubiff 



2 Gennady P. Berman, Alexander I. Nesterov and Alan R. Bishop 

2. The non-Hermitian approach to superconducting phase qubit 
measurement 

Non-Hermitian Hamiltonians naturally appear when the energy spectrum of a quan- 
tum system can be formally represented by both quasi-discrete (intrinsic) and con- 
tinuous parts, and one performs a projection of th e total wave function onto the 
quasi-discrete part of the spectrum l °M i| J l li l 1 ^ l 1 ' : > l j n this case , the corresponding 
intrinsic energy levels acquire finite widths which are associated with transitions 
from the intrinsic states to the continuum. Then, the dynamics of the intrinsic 
states can be descri bed by the Schrodinger equation with an effective non-Hermitian 
Hamiltonian CEHSl. 

We consider the current-biased qubit-continuum detector system based on a 
Josephson junction, which can be described by the Hamiltonian (h = 1) ^ : 



H t =^u> n \n)(n\+J2Prnn(t)\m)(n\+ f E\E)(E\dE +^ j a n (E)\E)(n\dE 

n myin n 

+ I a *n( E )\ n )( E \ dE (m,n = 0,1), (1) 

n J 

where w„ is the energy of level \n). The set of states {|n), \E}} forms a complete 
Hilbert space basis, so that an arbitrary state vector |?/>) can be expressed as 

|V) = C » + / C E \E)dE. (2) 

We use the following notation for qubit states: |0) = | \) and |1) = | f). The 
non-vanishing matrix elements, /3 TOn (i), are given by /3i2 = /?2i = cos(iot + tp), 
where f2o is the on-resonance Rabi frequency and w is the frequency of the biased 
ac current, I ac (t)- In general, the tunneling amplitudes, a n (E), are complex and 
depend on the total energy. 

Applying the Feshbach projection method 1 7 1 9 1 10 1 11 ^ we obtain an effective non- 
Hermitian Hamiltonian, H e ff = H — iW, which describes the dynamics of the qubit, 
where 

\ Slo cos(ut + ip) lo y ' 



and 



2 VTio To 

Here we consider the case in which the a„'s are real constants. Then, in terms of 
tunneling amplitudes, the rates of tunneling into the continuum are: Tq — 2-7rag, 
Ti = 2nal and Tio = r i = vToTi- 

The state of the qubit is controlled by a current I(t) = Id c — I^w cos(w< + ip) 
with a dc bias cur rent, Idc, and a time-varying bias current, I ac — I^w cos(uit + ip), 
with frequency, oj I ^U^ I The bias current, Idc, can be chosen so that the tunneling 
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is mostly from the upper level |1)SI. (See Fig. UJ) In what follows, we study the 
dynamics of quantum measurement by tunneling of a superconducting phase qubit 
in the limit Ti 3> To, 710, 7^, where 7^ is the dephasing rate. (We will discuss 
interaction with environment elsewhere lj.) 
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Phase difference 5 



Fig. 1. A superconducting phase qubit radiated by a microwave field; S is the gauge-invariant 
phase difference of superconducting order parameter across the junction; 710 denotes the rate of 
energy relaxation from |1) to |0) and Ti(ro) denotes the tunneling rate from level |1)(|0)). 



2.1. Rotating wave approximation 

In the rotating wave approximation, by removing fast oscillations, we find that 
Tieff fi, where 

- 1/Ao + A-iTi floe"* \ 

2^ n e^ X a -A-iT Q J- [) 

Here X = ujq + uji] A = lj w — oj is the detuning (with ui being the microwave 
frequency; and we denote wio = 0J\ — uiq). We neglect here the off-diagonal term 
r i, which describe the interaction between the two tunneling channels. We let 
the wave function \u(t)) = C\(t)\l) + Cb(i)|0) describe the qubit state. We assume 
below that (u(0)|m(0)) = 1. Then the solution of the Schrodinger equation with the 



a In a more complete model If, we consider the dynamics of the system within the generalized 
master equation in the Lindblad form, p = — i[H, p] + Lp— {W, p}, where {W, p} = Wp + pW and 
Lp = ^ ^)^_ 1 {2Lfe pL^, — L^Lkp — pL^Lk}, and transition operators, Lj,, describe the coupling 
to the environment 
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non-Hermtian Hamiltonian ([5]) can be written as 

fit „ . flt^ 



Ci (t) = e -^°'/ 2 U cos y - i cos 6 sin y ) d (0) - wT^ sin sin yC o (0) ) , (6) 
C7 (t) = e~ fXot/2 f ( cos y + t cos sin y ) C o (0) - ie 4V sin sin y Ci (0)") , (7) 



PiiW = e- r< 
Poo(t) - e _rt 



fit 

Y 

Of 



fit 
T 



cos — zcostf sin— ^)Ci(O) - ie ^sinflsin yC (0) 



fit 

fit 



/Sit iZt\ JZr 

(cosy + icos6>sinyjc (0) - ie^ sin (9 sin y Ci(0) 



(8) 
(9) 



Ti)/2, cos 61 



(A-i(r-r ))/O and sin (9 



where Ao = ujo + uj\ — iT, T = (Tq 

flo/fl, the complex Rabi frequency is 57 — y/fl^ + (A — i(T — To)) 2 ; the upper level 
population is pn(i) = |Ci(f)| 2 and the lower level population is poo = |Co(i)| 2 . 

In what follows we compare our theoretical results with experiments on quantum 
state tomography of a superconducting phase qubit and with the data obtained 
in the experiments observing coherent temporal oscillations between the quantum 
states of a Josephson tunnel junction In both experiments the tunneling rate 
from level |1) is more than 10 3 times that from level |0). So the total tunneling 
rate to the continuum is determined by the population of the upper level. The 
bias current was chosen so that the tun neling from level |0) is "frozen out" and 
escapes mostly from the upper level Therefore, the time-dependent tunneling 
probability is defined by the upper level population pn(t). 

In Fig. [2] we present our analytical solution for evolution of the Bloch vector, 
n(f) = (u(t)\cr\u(t)) . The decay rate, T, and the on-resonance Rabi frequency, flo, 
are taken from 3 : T = 0.035 ns" 1 , flo/{2n) = 80 MHz. Note, that the following 
relation between the density matrix and the Bloch vector, p — (nl -I- n • <x)/2, takes 
place. Our theoretical predictions shown in Fig. [5] are in a good agreement with the 
experimental data and with the numerical calculations reported i 

In Fig. [3] we compare our theoretical results with the reported experimental 
data ^ on evidence for the Rabi oscillations in a current-biased qubit based on the 
Josephson junction. The Rabi frequency, |f2|, and decay parameter, T, are chosen as 
inl = 8.9 • 10 6 rad/s and T = 0.204 /is" 1 . The initial state |u(0)) = 0.956|0) + 
0.291 1 1) , the detuning A, and the on-resonance Rabi frequency, ^o, are obtained by 
the fitting data to Eq. (|8]). We see good agreement between the experimental data 
and the predicted theoretical behavior. We find that our time-dependent upper- 
level population is a better fit to the experimental data than the theoretical results 
presented in Ref. ^by Eq.(5). However, additional information on the experimental 
data is required for better comparison. 

In the experiment reported in El, the total tunneling probability was identified 
as the upper level population given by, 



Pu{t) 



a 2 

M 2 



nt 



(10) 
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Fig. 3. The probability, P(t), of tunneling from the upper level to the continuum as a function 
of time. (n /(27r) = 0.47MHz, A/(2vr) = 1.34MHz, T = 0.204^s-\ To = 0.4- lCT 3 /us" 1 . The 
diamonds (red) are the data and the black solid line is the best fit to the oscillatory part of Eq. 



This is a particular case of our expression ([SJ , corresponding to the choice C\ (0) = 0. 



2.2. Tunneling into the continuum at zero ac current 

The fast frequency detector described in is based on a superconducting phase 
qubit tunneling into continuum. In particular, the measurement is performed with 
ac current, I^ w , equal to zero. Suppose, that initially the superconducting phase 
qubit is prepared as a superposition of ground and excited states. Writing \u{t)) = 
C\{t)\l) + Co(t)|0), we obtain the time-dependent amplitudes Co(t) and C\{t) as 
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solution of the Schrodinger equation 

<*<*) = e^((„« _ ™- si 4)a l( o) - ^ <*«,)), 

(11) 

C7 («) = e 0/ ( v (cos T + sm _jc o (0)-— sm T C 1 (0)j, 

(12) 



where il = y cj 2 — 2iu;io(r — To) — T 2 is the complex Rabi frequency. In the limit 
u>io 3> Ti, we have 



\u(t)) = e- l " ot e- r °* /2 Co(0)|0) +e-^ li e- rit/2 Ci(0)|l). 

■c perimen t 
f l6ll8ll9l 



(13) 



This simple formula is consistent with the known ex perime ntal data and the theo- 
retical description of a partial-collapse measurement ' 

The probability of escape from either |1) or |0) is given by 



Pesc{t) = l-Pll{t)- P 00(t), 

and in the same limit as above, ujiq 3> Ti, we obtain 

P esc {t) = 1 - /9 ii(0)e- ri * - (1 - pn(0))e- r °*. 



(14) 



(15) 



As one can see, the tunneling into the continuum is described by a double- 
exponential decay function P esc . This behavior of the probability of escape was 
observed in the experiment on measurement of dissipation-induced decoherence in 
a Josephson junction ^Sl. 

While the contribution of the off-diagonal term, Toi, in the experiments dis- 
cussed above is negligible, it can be more significant in the fast readout measure- 
ment of the phase qubit with a duration of measurement of approximately ~ 3ns 
l^JJ To estimate the contribution of the non-diagonal term, Tqi, we introduce the 
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Fig. 4. Deviation of population, F(t), as function of time. Left: Co(0) = 0, Ci(0) = 1. Right: 
C (0) = Ci(0) = 1/V2. 



function, F(t), which characterizes the difference in probabilities when Tqi ^ and 
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r M = 0: 



F(t) 



Pn(Toi,t) - pu(0,t) 



(16) 



In Fig.Uwe presen t th e results of computation for fast readout with the parameters 
being chosen as Ti = 0.1GHz, uj w /2t: = 5GHz and ri/r = 150. One can 

see that, depending on the initial conditions, the contribution from the interaction 
between the two tunneling channels can be up to « 0.1%. 

3. Concluding remarks 

In conclusion, we have presented a non-Hermitian description of the dynamics of 
measurement of a supeconducting phase qubit. We have shown that our theoretical 
predictions are in a good agreement with the experimental data obtained in various 
experiments on superconducting phase qubits. We believe that the approach based 
on non-Hermitian Hamiltonians can simplify the theoretical description of the dy- 
namics of the phase qubit measurement by tunneling to the continuum (in many 
cases and especially for a multi-qubit register). 

This work was carried out under the auspices of the National Nuclear Security 
Administration of the U.S. Department of Energy at Los Alamos National Labora- 
tory under Contract No. DE-AC52-06NA25396. The work by AIN was supported 
by research grant SEP-PROMEP 103.5/04/1911. The work by GPB was partly 
supported by the IARPA. 
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